LR and L+R Systems 



Bozidar Jovanovic 
February 10, 2009 



Abstract 

We consider coupled nonholonomic LR systems on the product of Lie groups. As 
examples, we study n-dimensional variants of the spherical support system and the 
rubber Chaplygin sphere. For a special choice of the inertia operator, it is proved that 
the rubber Chaplygin sphere, after reduction and a time reparametrization becomes 
an integrable Hamiltonian system on the (n — l)-dimensional sphere. Also, we showed 
that an arbitrary L+R system introduced by Fedorov in 15] can be seen as a reduced 
system of an appropriate coupled LR system Q 

1 Introduction 

In this paper we study nonholonomic geodesic flows on direct product of Lie groups 
with specially chosen right- invariant constraints and left-invariant metrics. 

Let Q be a n-dimensional Riemannian manifold Q with a nondegenerate metric 
•) and let I? be a nonintegrable (n — fc)-dimensional distribution on the tangent 
bundle TQ. A smooth path q(t) 6 Q, t 6 A is called admissible (or allowed by 
constraints) if the velocity q(t) belongs to T> q i t \ for all t € A. Let q — (qi, . . . , q„) be 
some local coordinates on Q in which the constraints are written in the form 

n 

(c4,?)=5>i*= ' j = l,-..,k, (1.1) 

i=l 

where a 3 are independent 1-forms. The admissible path q(t) is called nonholonomic 
geodesic if it is satisfies the Lagrange-d'Alambert equations 

k 

d dL dL v-^ i/ s . „\ 

^ = ^ + X>X(^, . = !,..,», (1-2) 

where the Lagrange multipliers Xj are chosen such that the solutions q(t) satisfy 
constraints (|1.1[) and the Lagrangian is given by the kinetic energy L = ^n(q,q) — 
\ Eij K ij<iiQj- After the Legendre transformation pi — dL/q\ = J]V KijQj, i = 1, ■ ■ ■ ,n, 
one can also write the Lagrange-d'Alambert equations as a first-order system on the 
cotangent bundle T*Q. As for the Hamiltonian systems, the Lagrangian L(q,q) (or 
the Hamiltonian H(q,p) — | Ysij ^PiPi m t ne cotangent representation of the flow) 
is always the first integral of the system. 

Suppose that a Lie group K acts by isometries on (Q, k) (the Lagrangian L is K- 
invariant) and let £q be the vector field on Q associated to the action of one-parameter 
subgroup exp(t£), £ £ 6 = TmK. The following version of the Noether theorem holds 
(see [U E]): if £q is a section of the distribution T> then 

Iff. *)-»■ M 
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On the other side, let £q be transversal to T>, for all £ € 6. In addition, suppose 
that Q has a principal bundle structure n : Q — > Q/K and that © is the collection 
of horizontal spaces of a principal connection. Then the nonholonomic geodesic flow 
defined by (Q,k,T>) is called a K-Chaplygin system. The system (|1.2[1 is Tf-invariant 
and reduces to the tangent bundle T(Q/K) = V/K (for the details see [26l [2l [8l [TT|V 

The equations (|1.2p are not Hamiltonian. However, in some cases they have a 
rather strong property - an invariant measure (e.g, see [T] 1271 [4]). Within the class 
of if-Chaplygin systems, the existence of an invariant measure is closely related with 
their reduction to a Hamiltonian form after an appropriate time rescaling dr = Mdt 
(see [lOl [291 [HI HE]). 

Veselov and Veselova [301 131] constructed nonholonomic systems on unimodular 
Lie groups with right-invariant nonintegrable constraints and left-invariant metrics, 
so called LR systems, and showed that they always possess an invariant measure. 
Similar integrable nonholonomic problems on Lie groups, with left and right invariant 
constraints, are studied in [171 1211 1221 O [19]. Recently, a nontrivial example of a 
nonholonomic LR system, which can be regarded also as a generalized Chaplygin system 
(n- dimensional Veselova rigid body problem [301 117] ) such that Chaplygin reducibility 
theorem is applicable for any dimension is given by Fedorov and Jovanovic [18] . 

It appears that LR systems can be viewed as a limit case of certain artificial systems 
(L+R systems) on the same group, which also possess an invariant measure (see Fe- 
dorov 15 ). The latter systems do not have a straightforward mechanical or geometric 
interpretation and arise as a "distortion" of a geodesic flow on G whose kinetic energy 
is given by a sum of a left- and right-invariant metrics. 

A class of L+R systems on G can be seen as a reduction of a class of nonholonomic 
systems defined on the semi-direct product of the group G and a vector space V (see 
Theorems 3, 4 in Schneider [28]). We shall prove that an arbitrary L+R system on G 
can be obtained as a reduction of a coupled noholonomic LR system defined on the 
direct product G x G. 

One of the best known examples of integrable nonholonomic systems with an in- 
variant measure is the celebrated Chaplygin sphere which describes a dynamically 
non-symmetric ball rolling without sliding on a horizontal plane and the center of the 
mass is assumed to be at the geometric center [9]. It is interesting that the Chaplygin's 
sphere appears within both constructions. In the construction described in [28] one 
should take for the configuration space the Lie group of Euclidean motion SE(3), that 
is the semi-direct product of SO(3) and R 3 [28] . On the other side, Chaplygin sphere 
is a LR system on the direct product 5*0(3) x R 3 (e.g., see [3]). This was a starting 
point in considering the coupled nonholonomic LR systems below. 

Outline and results of the paper. In Section 2 we recall the definition and 
basic properties of LR and L+R systems. We define the coupled LR systems and 
show that any L+R system can be obtain as a reduction of an appropriate coupled 
LR system (Sections 3, 4). An example of a coupled LR system on G x g is given, 
which provides an alternative generalization of the Chaplygin sphere problem (Section 
4, system (|6.17[) in Section 6). 

In Section 5 we study a n-dimensional variant of the spherical support system 
introduced by Fedorov [13] : the motion of a dynamically nonsymmetric ball S with the 
unit radius around its fixed center that touches N arbitrary dynamically symmetric 
balls whose centers are also fixed, and there is no sliding at the contacts points. 

Recall that the rubber rolling of the sphere S 2 over some other fixed convex surface 
in R 3 means that that in the addition to the constraint given by the condition that the 
velocity of the contact point is equal to zero, we have no-twist condition that rotations 
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about the normal to the surface are forbidden. The rubber rolling of the dynamically 
non-symmetric sphere over the another sphere, considered as a Chaplygin system on 
the bundle SO(3) x S 2 — > S 2 (where 5*0(3) acts diagonally on the total space), as well 
as the Hamiltonization in sphero-conical variables of S 2 is given by Koiller and Ehlers 
[12] . The integrable cases are found by Borisov and Mamaev [Jj. In particular, when 
the radius of the fixed sphere tends to infinity, we get the rubber rolling of the sphere 
over the plane (rubber Chaplygin sphere). The Chaplygin reducing multiplier for the 
rubber Chaplygin sphere is given in [11] . 

By the analogue, we define the n-dimensional rubber spherical support system with 
additional no-twist conditions at the contact points. It appears that both systems fits 
into the construction of coupled LR systems. Similarly as for the 3-dimensional spher- 
ical support system studied in [T3], we prove that the 3-dimensional rubber spherical 
support system is integrable (Section 5). 

Finally in Section 6 we consider the n-dimensional rubber Chaplygin sphere problem 
describing the rolling without slipping and twisting of an n-dimensional ball on an 
(n — l)-dimensional hyperplane TL in R n as coupled LR systems on the direct product 
SO(n) x R n_1 . It appears that the rubber Chaplygin sphere is a SO(n - 1) x R n_1 - 
Chaplygin system closely related to the n-dimensional nonholonomic Veselova problem, 
which allows as to prove the existence of the Chaplygin multiplier for a specially chosen 
inertia operator of the ball. In particular, when n = 3, the multiplier exist for any 
inertia tensor of the ball, and reduces to the one obtained in [TT1[T2] . 

2 Preliminaries 

LR systems. LR system on a Lie group G is a nonholonomic geodesic flow of a left- 
invariant metric and right-invariant nonintegrable distribution T> C TG (see [301 131j ). 
Through the paper we suppose that all considered Lie groups G have bi-invariant 
Riemannian metrics, or equivalently AdG-invariant Euclidean scalar products (•, ■) on 
corresponding Lie algebras. In particular, Lie groups G are unimodular. 

Let g = TidG be the Lie algebra of G. In what follows we shall identify g and g* 
by means invariant scalar product (•, •), and TG and T*G by the bi-invariant metric. 
For clearness, we shall use the symbol u> for the elements in g and the symbol m for 
the elements in g* = g. 

The Lagrangian is defined by L(g,g) = ^{Iuj,uj), where us = g -1 ■ g is the angular- 
velocity in the moving frame. Here I : g — > g is a symmetric positive definite (with 
respect to (■, ■)) operator. The corresponding left-invariant metric will be denoted by 
(■, )/. The distribution T> is determined by its restriction to the Lie algebra and 
it is nonintegrable if and only if is not a subalgebra of 0. Let f) be the orthogonal 
complement of with respect to (•, •) and let ai, . . . , a& be a orthonormal base of f). 
Then the right- invariant constraints can be written as 



(n,w = <^,fi 9 ) = o, 



() ff = Ad s -i(()) = ff - 1 -[)- ffj 



or, equivalently, 

(an,w)=0, a t = Ad s -i (a»), i = l,...,k. 

Here f2 = Ad 9 (aj) = g ■ g~ x represents angular velocity in the space. 
Equations (|1.2|l in the left trivialization take the form 



(2.1) 




(2.2) 



9 = 



g ■ w, 



(2.3) 
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where m = dL/doj — Iui £ g* is the angular momentum in the body frame. 

The Lagrange multipliers A,, can be found by differentiating the constraints (J2HJ. 
They are actually defined on the whole phase space T*G and we can consider the 
system (|2.2[) . (|2,3|l on T*G as well (see [31]). The constraint functions {cti,u>) are 
then integrals of the extended system and the nonholonomic geodesic flow is just the 
restriction of (|2.2[) . f|2.3|) onto the invariant submanifold (|2,1[) , 

Instead of (|2.2[1 . (|2.3[) . one can consider the closed system consisting of f|2.2|) and 

ati = [ai,u], i = l,...,k, (2.4) 

on the direct product g 1+fc = a 1 , . . . , a fe )}. Let I^ 1 \na = pr^g o/ _1 o pr^ , where 
pr^ B is the orthogonal projection to f) 9 . Then the system (|2.2p . (12. 4p has an invariant 
measure with density n — \J det(/ _1 |f>s ) = ^/det ((J _1 (ai), ctj)) (see [31]). 

Also, since for ( 6 j, the associate vector field £g of the left G-action is right in- 
variant and the momentum mapping of the left action equals to M — Ad 9 (m) (angular 
momentum in the space), the LR system (|2.2|l . (|2.3[1 has the Noether conservation laws: 

^(Ad 9 (m),O=0, (2.5) 

If the linear subspace f) is the Lie algebra of a subgroup H C G, then the Lagrangian 
L and the right-invariant distribution T> are invariant with respect to the left //-action. 
As a result, the LR system can naturally be regarded as a i!f-Chaplygin system |18j . 

Geodesic flow on G with L+R metric. In addition to the nondegenerate linear 
operator / defining the left-invariant metric introduce a constant symmetric 

linear operator IT : g — > g defining a right-invariant metric (•, -)n on the n- dimensional 
compact Lie group G: for any vectors u, v 6 T g G we put (u,v)n = (ug" 1 , T^vg^ 1 ). 
We take the sum of both metrics and consider the corresponding geodesic flow on G 
described by the Lagrangian 

L = ^>-M + Tjig^g' 1 ,n° ' gug' 1 ) = -(w,Iu>) + {w,n 9 uj}, 

where IF = Ad g -in Ad g . We can also consider the case when II 9 is not positive 
definite, but the total inertia operator B = I + II 9 is nondegenerate and positive 
definite on the whole group G. 

The geodesic motion on the group is described by the Euler-Poincare equations 

r i -idL dL . . 

m=[m,u\ + g —-, m = — =Buj, (2.6) 
og ouj 

together with the kinematic equation g = g - w. 

In order to find explicit expression for g~ 1 (dL/dg), we first note that for any £ € g, 
(£, g~ 1 (dL/dg)) — v^(L), where is the left-invariant vector field on G generated by 
£. Since the metric (•,•)/ is left-invariant, we have 

1 1 T 

v^(L) — — v^((lj, Huj)) — —{lo, iladjw + ad j Hlj) = (TIlj, [£, cJ\) = (£, ad w Hit)}. 

As a result, g~ 1 (dL/dg) = ad u Ilw. 

Also, in view of the definition of II, its evolution is given bynxn matrix equation 

II = IIad w 4-ad^IL (2.7) 

Since (•, ■) is AdG invariant scalar product, we have ad^ = — ad„, and II = [II, ad w ]. 

Equations (|2,6[1 , (|2.7[l form a closed system on the space g x Symm(n) with the co- 
ordinates u)i,Hij (u> — £\ ^i e i, n = Yli<j Tlij&i'Si&j), where e\, . . . , e n is a orthonormal 
base of g. 
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L+R, systems. Following Fedorov 15 , consider the equations (|2.6p modified by 
rejecting the term g~ 1 (dl/dg). As a result, we obtain the another system 

[Bw,w], g = g-u>, B = I + U (2.8) 

[Bu,u], ^n = nad a; +ad^n (2.9) 

on the space g x Symm(n) . This is generally not a Lagrangian system, and, in contrast 
to equations (|2.6p . (|2.7|l . it possesses the "momentum" integral {Buj,Buj}. In view of 
the structure of the kinetic energy, we shall refer to the system (12. 8[) (or (|2.9|0 as L+R 
system on G [15] . 

The L+R system (|2.9|l possesses also the kinetic energy integral h{uj,Bu)) and an 
invariant measure (in coordinates Wj, Uy) pdidi A ■ ■ ■ A a; n A dlln A ••• A II nn with 
density ^ = ^/det(/ + n) (see [T41IT5] ). 

As mentioned above, a nonholonomic LR system on a Lie group G can be obtained 
as a limit case of a certain L+R system on this group. Indeed, suppose that the 
operator defining a right-invariant metric on G is degenerate and has the form LT = 
e(ai ® ai + ■ ■ ■ + at, ® o^), k < n, e — const > 0, where, as in (|2.1[) . ai, . . . , a*, are 
orthonormal right-invariant vector fields cu = g" 1 ■ en • g, cii = const £ @. The L+R 
system (|2.9p on the space (w, ai, . . . , a*,) can be represented in form 

Jci = /(Z + n)" 1 ^,^, II = IIad u +ad^n. (2.10) 

Then the following statement holds (see |15|). As e — > oo, the equations (|2.10p 
transform to the equations with multipliers (|2.2p and constraints (|2.ip . where m = to. 



on TG, or the system 



Tt = 



3 Coupled nonholonomic LR Systems 

Define a coupled nonholonomic LR system on the direct product Gxft (G = Gi) as 
a LR system given by the Lagrangian function 

L = i(to,o;> + i J D(w,w) (3.1) 

and right-invariant constraints 

(n,& )=0, (3.2) 
(fi + Pi W, f)j) = 0, i = l,...,g, (3.3) 

where f)i, i = 1, . . . ,q are mutually orthogonal linear subspaces of q. 

Here (w,w) = (<? _1 <?,S 1 ~ ffi), is the angular velocity in the body and (fi, W) = 
Ad( 9iS1 ) = (Adg(oj), Ad 91 (w)) is the angular velocity in the space, (g,gi) 6 G x Gi. 
The constant D is greater than zero, while pi, i = 1, . . . ,q are arbitrary non-zero, real 
parameters. 

The Lagrangian (|3.ip in the second variable is right-invariant as well. It is conve- 
nient to write the equations of motion both in the left-trivialization (in variables g and 
uj) and right-trivialization (in variables g\ and W) 

T(G x Gi) « G x Gi x x 01 = {{g, 9i ,lo,W)}. (3.4) 

Then the right-invariant distribution T> C T(G x ft) is given by 

P = {(<?, Pl , ^, W) | (Ad 9 H,f,o)=0, (Ad s (^) + Pl W,t )l > =0, i = l,...,g}. 

Let f)f = Ad g -i (f)i) = g~ x -i)i-g and let pr^s : g — » f)f be the orthogonal projections, 
i = 0, 
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Proposition 3.1. The admissible path (g{t), gi{t), u)(t), W(t)) is a motion of the non- 
holonomic LR system \'J. fS.HAl , \3. S[j if it satisfies equations 

-^rVsg^aMu,]), (3.5) 

(Ad 9 u), (3.6) 



Bu 


— [Ilo,lo] - 


DW 


--±7 

i=l Pl 


9 


= g ■ w, 







(3.7) 
(3.8) 

tu/iere S = J + n = J + J^Li ^V/ 9 ? P r ij? and = P r ijg oS_1 ° P r i,g : -» fto- 

Proof. The equations of a motion in the right-trivialization (or in the space frame) 
read 



M = 


i=0 


(3.9) 


DW = 


<3 
i=l 


(3.10) 


g = 


Q ■ g, 


(3.11) 


gi = 


W-g u 


(3.12) 



where the Lagrange multipliers (reaction forces) Aj belong to fjj (i = 0, 1, . . . , g) and 
M = Ad 9 (/oj) is the first component of angular momentum in the space frame (the 
second component is M\ — DW). 

Differentiating the constraints (13.3|) . from (|3.10p we obtain 



3=1 

that is 

A* = 2-pr & .(f2), i = l,...,g. (3.13) 

The equation (|3.6[1 follows from (|3.10|l . (|3.13|l and the relation 

Q, = Ad 9 Cj. (3.14) 

From (|3.13[l and identities (|3.14|l . pr^s = Ad 9 -i pr h . Ad 9 and 

Ad s -i M = Id) + [u, Iuj], 

the equation (|3.9|l in the left-trivialization takes the form 

Buj = [7w,w] + A , A = Ad g -i(A ). (3.15) 

Now it remains to find the Lagrange multiplier Ao- Differentiating (|3.2|l we get 

(ft, ft) = (Ad 9 (w),f)> = <w,bg)=0. 

Whence, according (|3.15|l it follows Ao = — (S -1 | fc s) -1 pr^s S~ 1 ([/tj,tj]). The proof is 
complete. □ 

The Lagrangian (|3.1[) as well as constraints (|3.3|l are right ({Id} x Gi)-invariant 
and the equations (|3.5|l . (|3.6|) . (|3.7|l can be seen as a reduction of the system to 



V = ©/({Id}xG0 = {( S ,w,W) | (Ad 9 (o;),fto> = 0, (Ad B (w)+ Pi W, ^) = 0, t = 1, ...,«}. 
Let T>o C TG be the right-invariant distribution defined by (|3.2|l . 
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Theorem 3.2. The equations fS. 5\) . f3. 61) . \3. 7| ) on T> are reducing to the following 
system on T>o C TG: 

j t {Bu) = [BuM-{B-Xir 1 VhlB-\\I"M\ 9 = 9-u>. (3.16) 

Proof. The equations (|3.5[) and (|3.7|l form a closed system on T>q. If (g(t),ui(t)) is a 
solution of (13. 5[) . ()3.7|) . then one can easily reconstruct the motion of W. Let 



6= (f)i + ■■■ + [),)• (3-17) 

From l|3.7j> we have 

^pr e W = 0, (3.18) 

while the (^-components of the angular velocity W are determined from the constraints 
0: 

P^, W = -1/PiPttn Ad 9 (tj), i = l,...,q. 
Now, let oi, . . . , afcj be the orthonormal base of hj . Then oti = Ad fl -i (ai), . . . , = 
Ad g -i(afe.) will be the orthonormal base of h?. We have 



£=1 i=l 

Whence, by using (12. 4p and the identity {lu, [cti,ui]} — 0, we obtain 



^pr(,9(w)j = ^ ((ij,Qj)ai + (u, [q,,cj])q, + (uj,ai)[ai,ui]) 



= pr 6 »(w) + ^(w,a i )[a i ,w]. 

i=l 

The above equation implies that (13. 5[) . (|3.7p can be rewritten in the form ()3. 16[) . □ 

The derivation of {Bui, lu) along the flow is: ^ {Bui, ui) = 2{[B, lu], uj) + 2{Ao, lu). The 
first term is equal to zero since (•, ■) is a Adc-invariant scalar product, while the second 
term is equal to zero from the constraint (|3.2[) . We can refer to L re£ j = ~{Bui, u>) as to 
the reduced Lagrangian, or reduced kinetic energy. If pr E W = 0, the reduced kinetic 
energy coincides with the kinetic energy of the reconstructed motion on the whole 
phase space. 

From the equation (|3.9p we also get the linear conservation law 



j (pr (o Ad fl (Iw)) =0, where { = (f)o + f)i + • • • + f)?)^ • (3.19) 

The integrals (|3.18[) and (|3.19[) are actually Noether integrals (|2.5|l of the system. 
The other Noethers integrals are trivial: 

fpr,,. Ad 9 (Jw) - j- pr b . W^j = 0, i = l,...,q. 

Remark 3.1. If f)o = 0, i.e., we do not impose the constraint (|3.2p . the reduced system 
is an L+R system on the Lie group G 

l(Bu) = [Bu,u], 9 = 9-u. (3-20) 
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Further suppose that (|3.17|) is the Lie algebra of the closed Lie subgroup K C G and 
that linear subspaces hi are Ad^ -invariant: 

Ad fc hi = hi, k e K, i = l,...,q. 

Then, since f)* 9 = hf, k G K, the L+R equations (|3.20p are left Tf-invariant and we 
can reduce them to Q x q, where Q — G/K is the homogeneous space, with respect to 
the left-action of K. 

Remark 3.2. In the case when h is the Lie algebra of a closed subgroup H C G, 
f)i + f)2 + -- - + f)ij = and linear spaces f)j are Adjj invariant, then the coupled LR 
system (13. (|3.2|l . (|3.3|) is (// x Gi)-Chaplygin system with respect to the action: 

{a, b) ■ (g,gi) = {ag^^' 1 ), (a,b) G H x Gi . 

The reduced space £>/ (_ff xGi) is the tangent bundle of the homogeneous space G/H. 

Theorem 3.3. An arbitrary L+R system \Z. 80 can be seen as a reduction of an ap- 
propriate coupled LR system. 

Proof. Let e\ , . . . , e n be the orthonormal base of q in which the symmetric operator 
n° has the diagonal form: II = X^™=i ai e i® e i- Then the right invariant term in (|2.8[1 
reads II = II 9 = Ys?=i Gi e f ® e f ; wnere e f are given by 

e? = Ad B -i (ei), . . . , e 9 = Ad s -! (e n ). (3.21) 

Consider the 2-coupled nonholonomic LR system (|3.1[) , (|3.3|) . where q = n and f)i 
are the lines in the directions of ej, i = 1, . . . , n. We can choose parameters D, pi, such 
that Oi — D/pf , i = 1, . . . , n. The system represents a {Id} x Gi- Chaplygin system 
with reduced equations of the required form (|2,8|l , □ 

4 A^-Coupled Systems 

There is a straightforward generalization of the construction to the case when we have 
coupling with TV different Lie groups, that is the configuration space is the direct 
product G x Gi x • • ■ x Gjv and the Lagrangian is 

1 1 N 

L= -{Iw,w) + -Y,Di(wi,v/i)i, (4.1) 

~ i — 1 

where (•, -)j are Ado ; invariant scalar products on Lie algebras @; = T^G;, i = 1, . . . , N. 
Let us fix a base e\, . . . ,e n of g and some bases fx,. . . , fd t of $ji (d; = dimg;). Let 

Ai :g^R Pi , B t : Si -» » p % i = l,...,JV. 

be the linear mappings with matrixes [j4»] (pi x n) and [Bi] (pi x di) in the above bases. 
In addition, we suppose that the (pi x pi)-matrixes 

[d] = [Bi][Bi] T , i = l,...,N 

are invertible. Consider the right invariant constraints given by 

A i Q + B i W i = 0, i = l,...,N. (4.2) 

Here, u, w, and Q,, Wj are velocities in the left and right trivializations, respectively 
and Di > 0, i = 1, . . . , N are real parameters. 
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Let [Q], [WJ denote the column matrix, representing fl and W.; in the chosen 
bases. We have [uj] g = [Q], where [£] 9 is the column, representing £ 6 g in the base 

In the right-trivialization, the equation in Wj reads 

Di[Wi] = [Bi] T [\i], (4.3) 

where [Aj] is the Lagrange multiplier (pi x l)-matrix. Differentiating the constraints 
P~2")) . from (JO) we get 

[X i ] = -D[Ci]- 1 [A i ][a], i = l,...,N. (4.4) 

Repeating the arguments of Theorems 13.11 and 13.21 the considered Af-coupled non- 
holonomic system reduces to the L+R system 

^-(Bu) = [Bw,<j], g = g-uj, 

where Blo = Iui + Ilea, and Ulo in the matrix form, relative to the base (|3.2ip . is given 
by 

N 

^u] g = ^2Di[A^ T [Ci]- 1 [A i ][u) B . 

i=l 

As above, one can easily incorporate an additional right invariant constraint of the 
form ([32}. 

LR systems on G x g x • • • x g. As an example, consider the case where d are 
all equal to the Lie algebra q considered as a Abelian group, (■, = (■, ■) and the 
constraints (|4.2|l are given by 

[n,r i ] + p<w i = o, i = i,...,N, (4.5) 

where Ti are fixed elements of the Lie algebra g and pi are real parameters. Note that, 
since Gi = g is Abelian group, the angular velocities coincide with the usual velocity: 
£i = W, = Wj, £ e g. 

The equations of a motion in the right-trivialization read 

JV 

M = ^2[JU,T t ], g = V-g (4.6) 

j=0 

D,W, = Pi Ai, £i = W h (4.7) 

where M = Ad g {Iu>). This is a {Id} x g^-Chaplygin system and it is reducible to TG. 
Differentiating the constraints (|4.5|) . from (|4.7|l we get the Lagrange multipliers 

Ai = -^[Ti,n], i = l,...,N. 

Pi 

Therefore, the equations (|4.7[) in the left-trivialization take the form 

N Di 

7w=[Jw,w]-^-j[[7i,w],7i], g = g-uj. 

i=l Pi 

where ji = Ad g -i (Ti), i — 1, . . . , N. Next, from the identities 

^[[7*5^,7^ = [[7<,w],7i] - [[[-yi,u;],ji],Lj], i = l,...,N, 
we obtain the following proposition 
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Proposition 4.1. The reduced equations of the N -coupled nonholonomic system {4.1}) , 
(4-5p are given by the L+R system 

j t {Bu) = [BuM, 9 = 9-", (4-8) 

where 

Bu) = IW + -j[[7iiW],7i]- 

Remark 4.1. Nonholonomic systems on semi-direct products G x CT V, where a is a 
representation of the Lie group G on the vector space V are studied in Schneider [28] . 
Proposition 4.1 can be derived from Theorem 3 given in |28j . 



5 Spherical Support 

Consider the motion of a dynamically nonsymmetric ball S in R n with the unit radius 
around its fixed center. Suppose that the ball touches N arbitrary dynamically sym- 
metric balls whose centers are also fixed, and there is no sliding at the contacts points. 
We call this mechanical construction the spherical support. For n = 3 spherical support 
is defined by Fedorov [13 1 115] . 

The configuration space is SO(n) N+1 : the matrixes € SO(n) map the frames 
attached to the ball S and the ith ith peripheral ball to the fixed frame, respectively. 
The Lagrangian is of the form Q4.ip . where for (•, •) we take the scalar product propor- 
tional to the Killing form 

(X,Y) = -~tr(XY), (5.1) 

the angular velocities u), fi, w,, Wj of the balls are defined as above, / : so(n) — > so(n) 
is the inertia tensor of the ball S and Di,pi g R are the central inertia moment and 
the radius of the ith peripheral ball. 

Let Fi £ R n be the unit vector fixed in the space and directed from the center C of 
the ball S to the point of contact with the ith ball. Nonholonomic constraints express 
the absence of sliding at the contact points. This means that velocity of the point of 
contact of the ball S with the ith ball, in the space frame, is the same as the velocity 
of the corresponding point on the ith ball. 

Consider the fixed point on the ball S with coordinates r and R in the body and 
space frames, respectively. Then the velocity of the point r in space is given by the 
Poisson equation (e.g, see [17]) V = R = -^(g • r) = g ■ g^ 1 ■ g ■ r = QR. Therefore, the 
velocity of the contact point with the ith peripheral ball is given by S7P; . Similarly, the 
velocity of the corresponding contact point of the ith ball in the space frame is given 
by — piWTi and the constraints are 

fir, + piWiTi = 0, i = l,...,N. (5.2) 

We see that the n-dimensional spherical support is actually a TV-coupled LR system 
studied in the previous section. Let 

l l =g~ 1 T l , i = l,...N (5.3) 

be the contact points of S with the ith ball (i — 1, . . . , N) in the frame attached to the 
ball S. Then the right-invariant constraints (|5.2|l can be rewritten in the form 

(n + Pi Wi,t)i) = 0, i = l,...,N, (5.4) 

where 

d 1 =R n Ar l , f)l = Ad 9 -i(b0 =» n A 7 *, i = l,...,N 
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Figure 5.1: The spherical support 



are linear (no mutually orthogonal) subspaces of the Lie algebra so(n). 

From the identity pr^s d> = (cL>"/i) A"/i — ill 74 ® 74 + 7, <8 74 oj, the equations of the 
motion become 

N A 

Iio — [Iuj, lu] - 22 —5 (w 7» ® 7» + 7i ® Ji w) , g = g-u 

AWi = -- fnr i ®r i + r j <8)r j n') , gi = w z - gi , % = i,...,n. 

Pi \ / 
We have the conservation laws 

w,i-Wiri(8ri-ri8riWi = o, i = i,...,jv, 

which together with the right ({Id} x S'0(n) iv )-symmetry lead to the following state- 
ment 

Proposition 5.1. TTie spherical support system reduces to the L+R flow 

^(Bw) = [Bw,w], g = g-u, (5.5) 

where Bw — Iuj + Y2i=i Di/Pi {u> 74 ® 74 + 74 <8> 74 w) and 74 are defined by f5.3J) . 

One can say that the reduced system (|5.5[) on TSO(n) describes the free rotation 
of a "generalized Euler top", whose tensor of inertia is a sum of two components: one 
is fixed in the body and the other one is fixed in the space. 

Note that the vectors 74 in the frame attached to the ball S satisfy the Poisson 
equations (e.g., see [17] ) 

74 = -u>7i, i = 1, . . . , N. (5.6) 
By introducing Xt = 74 ® 74, from (|5.6[) we obtain 

Xi = lXi,uj], i = l,...,N. (5.7) 
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Combining (|5,5|l and (|5.7p we get family of integrals - the coefficients of the poly- 
nomials 

JV 

tr{Bu + J2n i X i ) k , k = l,...,n. (5.8) 

i-l 

For n = 3 the system is integrable by the Euler-Jacobi theorem, and its generic 
invariant manifolds are two-dimensional tori (see [131 115| ). 

Remark 5.1. If the positions of peripheral balls are mutually orthogonal 

(r^Tj) = ( 74)7j ) = Sv, l<i,j<N< n, 

then the components of ji can be seen as redundant coordinates on the Stiefel variety 
V(n, N) = SO(n) /SO(n - N). The system is invariant with respect to the SO(n - N) 
action, representing the rotations in the space orthogonal to span{7i, . . . ,7jv}. The 
SO(n - 7V)-reduced system on TSO(n)/SO(n — N) = V(n,N) x so(n) is given by 
Poisson equations (|5,6|l and the first equation in (|5,5ll . 



Rubber spherical support. Now consider the rubber spherical support system 
in R n . The analogue of rubber rolling is that, in addition to the constraints (|5.4|) , the 
rotations of the ball S and ith peripheral ball around the vector Fi are the same: 

(n-W i ,l i ) = 0, i = l,..., AT, (5.9) 

where 

li = f)i, ti = so(n-l). 
Since pr e . = I — pr^ . we get 
Proposition 5.2. The rubber spherical support system is described by the equations 

^(B*w) = [B*w,w], (5.10) 

ftWj = Ditl - D^—^- [tlYi ® r t + r t ® Ti f2 j , i = l,...,N, (5.11) 

g = g-iu, (5.12) 
g i = W i -g i , i = l,...,N, (5.13) 

where 

N 1 - p 2 

!=1 Pi 

The equations (15.11[1 are trivial since W can be expressed in terms of £1 from 
constraints (|5.4|l and (|5.9[) . 

As above, we get family of geometric integrals that can be expressed as the coeffi- 
cients of the polynomials 

N 

tr{B*uj + ^^Xi) k , k = l,...,n. (5.14) 

i—l 

For n = 3, among the reduced kinetic energy i(B*w,oj) and integrals (|5.14[) there 
are four independent one. 

Theorem 5.3. Forn = 3, the rubber spherical support system (5.1 0)) . 15.1S\) is solvable 
by the Euler-Jacobi theorem and its generic invariant manifolds are two-dimensional 
tori. 
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6 Rubber Chaplygin Sphere 



Following [171 116j . consider the generalized Chaplygin sphere problem of an n-dime- 
nsional ball of radius p, rolling without slipping on an (n — l)-dimensional hyperspace 
H in R n . For the configuration space we take the direct product of Lie groups SO(n) 
and R n , where g £ SO(n) is the rotation matrix of the sphere (mapping frame attached 
to the body to the space frame) and r £ R n is the position vector of its center C (in 
the space frame). For a trajectory (g(t),r(t)) define angular velocities 

u = g~ 1 g, fl^gg' 1 , w = W = r. 

The Lagrangian of the system is then given by 

1 1 

L = -(Iu,^} + -m(w,w). (6.1) 

Here I : so(n) — > so(n) and m are the inertia tensor and mass of the ball, (•, ■) is given 
by (15. If) and (•, •) is the Euclidean scalar product. 

Let r £ R™ be a vertical unit vector (considered in the fixed frame) orthogonal to 
the hyperplane TL and directed from TL to the center C. The condition for the sphere 
to role without slipping leads that the velocity of the contact point is equal to zero: 

-pfir + W = 0. (6.2) 

This is a right-invariant nonholonomic constraint of the form (|4.2[) . If we take 
the fixed orthonormal base Ei = (1, 0, . . . , 0, 0) T , . . . , E n — (0, 0, . . . , 0, 1) T , such that 
r = E n , then the constraint (|6.2p takes the form 

h = p^lin, i = l,...,n — 1, r n = 0, where fijj = (fi, Ei A Ej). 

The last constraint is holonomic, and for the physical motion we take r n = p. From 
now on we take SO(n) x R n_1 for the configuration space of the rolling sphere, where 
R™~ is identified with the affine hyperplane pF + TL. 

Let f) C so{n) be the linear subspace f) = R n A F and I = so{n — 1) its orthogonal 
complement in so(n). Define the rubber Chaplygin sphere as a Chaplygin sphere (|6,1[1 . 
(|6.2p subjected to the additional right-invariant constraints 

(fi,t) = (w,6 9 ) = 0, e 9 = Ad s -ie, ^-=0, l<i<i<n-l, (6.3) 

describing the no-twist condition at the contact point. As a result, the distribution 

V = {(g,r,u,W) | (u,t s ) = 0, W = P Ad 9 Hr}. 

is right SO{n) x R n_1 as well as the left SO{n - 1) x R n_1 invariant (SO{n - 1) is the 
subgroup of SO(n) with the Lie algebra t). Moreover, the rubber Chaplygin sphere is 
a (SO(n - 1) x R n ~ 1 )-Chaplygin system. 

Let 7 be the vertical vector in the frame attached to the ball 7 = g~ l T. Then 

[j 9 = Ad s _i (t>) = I™ A 7 =: f| 7 , e fl = Ad g _i (6) = (R n A 7) X =: f 

and the reduced space V/(SO(n - 1) X R n_1 ) is the tangent bundle TS" 1 " 1 of the 
sphere which can be identified by the position of 7. 
The equations in the right-trivialization read 

M = A - pAi A T, g = f2 • g, (6.4) 
mW = Ai, r = W, (6.5) 
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where M = Ad g (Iui) is the ball angular momentum in the space and Ao £ f), Ai € R n 
are Lagrange multipliers. From (|6.2p and (|6.5|l we find Ai = mpClT. On the other hand 

Ai AT = mp(tir) AT = mpftir ®T + T <g>rn\ = mppi^(Cl). (6.6) 

Whence, we can write equations (|6.4|l as a closed system on Do C TSO(n), where 
T>o is the right-invariant distribution defined by (|6.3p (reduction of R n_1 -symmetry). 
From (|3.14|l . (|6.6|l and the relation pr^ (d>) = (cj • 7) A 7 = 7 ® 7 + 7 ® 70;, in the 
left-trivialization of TSO(n) the reduced system takes the form 

Ilu = [Jw,a;] — mp 2 (dj 7®7 + 7®7tj) + Ao, g = g ■ u>, 

where Ao = Ad fl -i(Ao). Let 

k = Ilu + mp 2 pr,j T w = I id + mp 2 (u) 7(87 + 787^) G so(n)* (6-7) 

be the angular momentum of the ball relative to the contact point (see [IT]). Then we 
have: 

Proposition 6.1. The motion of the rubber Chaplygin sphere, in variables Lu,g, is 
described by 

k=[k,^] + A , g = g-u, (6.8) 
or, in variables lu,j, by equations 

k=[k,w]+Ao, 7 = -W7. (6.9) 

The Lagrange multiplier matrix Ao belongs to t 1 and is determined from the constraint 
fOI) . 



Reduction and Hamiltonization. From the constraints (|6.3|) . the momentum 
(|6.7[) equals to k = Ilu + mp 2 to. Therefore, as in the 3-dimensional case [111 [6], the 
equations (|6.8|l are equivalent to the motion of a rigid body about the fixed point with 
the left-invariant kinetic energy given by the inertia operator I + mp 2 I and constraint 
(|6.3p (n-dimensional Veselova rigid body problem [3011 17] ). 

Now we simply follow [15]. The reduced Lagrange-d'Alambert equations of the 
rubber Chaplygin sphere (|6Tj) , lfO]l, (|Q)l on 



r5 ,n-l ^ 2y( SO ( n _ !) x R »-l) V /SO(n - 1) 

are given by 

=« / + V)*(7,7),(I-pr () ,)[*(7,7),*(7,0]>, (6-10) 

for all virtual displacements £ 6 T-yS" 1-1 (see [18]). Here $(7,7) = 7 A 7 is the 
momentum mapping of the right 5*0(n)-action on the round sphere S" l_1 and the 
reduced Lagrangian is given by 

Lr ed (l,i) = ^(U + ^P 2 )*(7,7),*(7,7))- (6.H) 
After the Legendre transformation 

dL red dL reg <9$ 2 . rA 

we can also write the reduced Lagrange-d'Alambert equations as a first-order system on 
the cotangent bundle t*^™- 1 which is realized as a subvariety of R 2n = (q,p) defined 
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by constraints (7,7) = 1, (7,p) = (since 7$ is skew-symmetric, the momentum p 
satisfies (7,p) = 0). The system takes the symmetric form 

7= _$( 7j7 ( 7j p)) 7i p - -$(7,7(7, p))p, (6.12) 

where 7 = 7(7, p) is the inverse of the Legendre transformation. 

Let a be the canonical volume 2(n — l)-form on T*S n ~ 1 . Then we have (see [18]): 

Proposition 6.2. The reduced system j6.12\) on T*S n ~ possesses an invariant mea- 
sure 

l/Vdet(J + rap 2 I|i,7) a, I + mp 2 I\t,-, = pr^ o(I + mp 2 I) o pr^ . 

Furthermore, as it follows from |18| . with the operator / defined on the bi- vectors 
X A Y by a diagonal matrix A = diag(j4i, . . . , A n ) by 

I(X A Y) = AX A AY-mp 2 X AY, (6.13) 

the Chaplygin reducibility is applicable for any dimension: 

Theorem 6.3. (i) If the inertia operator is given by S6.13\) . the density of an invari- 



ant measure in Provosition \6.2\ takes the following simple form 

(A 7 ,7)- ( - 2)/2 . 



(ii) Under the time substitution dr — 1/ ' \J (A7, 7) dt the reduced system f6.1(fy (or 
16.12]) ) becomes a Hamiltonian system describing a geodesic flow on S" I_1 with 
the Lagrangian 



7, 



(*SH[("&£W>-(*-S) 



(6.14) 



(Hi) For A with distinct eigenvalues, the latter system is algebraic completely integrable 
and generic invariant manifolds are (n — I) -dimensional tori. 

(iv) Moreover, the SO(n — 1) -reconstruction of the motion is solvable: the generic 
trajectories of the system 16. 8\ ) are straight-lines (but not uniform) over (n — 1)- 
dimensional invariant tori. 

The complete integration is presented in [18] . Given a solution (g(t),uj(t)) of the 
system (|6.8p , the reconstruction of r- variable simply follows from the integration of the 
constraint (16.211 



r(t) - r(t ) = p / Ad 9(t) u(t)T dt . 
In the case n — 3, under the isomorphism between so(3) and R 3 

LOij = EijlLJl, kij=£ijlkl, (6.15) 

from (|6.9[) we obtain the classical rubber Chaplygin's ball equations [TT] 

k = k x lo + A7, 7 = 7x0}, (6.16) 

where A is determined from the constraint (io, 7) = and k = ICj+mp 2 id—mp 2 (uj, 7)7 = 
Ilo + mp 2 uj. 

For n = 3, the relation (|6.13|l defines a generic inertia tensor. Thus the rubber 
Chaplygin sphere in R 3 is integrable. Indeed, let / : R 3 — ► R 3 be an arbitrary inertia 
tensor. Under the isomorphism (I6.15|l . the matrix A is determined from (|6.13|l via: 



A = A(J + mp 2 I)~ , A = N /det(7 + mp 2 I) 
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The Hamiltonizaton of the reduced system on T* S 2 is obtained in [111 112] . The 
Chaplygin multiplier given in Theorem 16,31 

dr = dt/^A((I + mp 2 I)- 1 1 ,~f) 

up to the multiplication by a constant, coincides with the expression obtained in [111 

□2]. 

Remarks on the Chaplygin sphere. 

• Note that the Chaplygin sphere equations 

k = [k, cj], 7 = — uj'y, k = Ilo + mp 2 (^718)7 + 7(8)70;) 

coincide with the reduced equations of the spherical support system for N — 1, 
where instead of Di/px we should put mp 2 . This is not the case for rubber 
analogues of the systems. 

• Borisov and Mamaev [5j [6] proved that the classical Chaplygin rolling sphere prob- 

lem is Hamiltonian after an appropriate time rescaling. Recently, the Hamil- 
tonization of the homogeneous Chaplygin rolling sphere problem in R n is given 
in [20] . while the Hamiltonization of the non- homogeneous reduced Chaplygin 
sphere problem is obtained in |25j . 

• Let us turn back to the coupled LR system described in Proposition 4.1. Take N = 1 

and denote Ti = T, 71 = 7, D\ = m, pi = 1/p. The system (|4.8[l is additionally 
Gr-invariant, where Gr C G is the isotropy group of T. Let O = G/Gr be the 
adjoint orbit of T. Then (|4.8|) reduces to O x g: 

k=[k,w], 7 = [7,w], k = Blo — Ilo + mp 2 ^, w],7]. (6-17) 

For G = SO(3) we reobtain the equations of a motion of the Chaplygin sphere 
in R 3 . Thus the system (I6.17P can be seen as an alternative generalization of the 
Chaplygin sphere problem. 
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